Reduction of superfiuid gap by scattering 
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The effect nucleon dressing by scattering on the pairing gap in nuclear matter is discussed. Numer- 
ical results from self-consistent T-matrix calculations are compared to quasi-particle approximations. 
The dominant effect of scattering can be accounted for by a renormalization of the quasi-particle 
strength. The exponential sensitivity of the pairing gap and the critical temperature on the pairing 
strength makes it strongly dependent on many-body effects. 
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I. INTRODUCTION 



' Nuclear systems at low temperatures undergo a superfiuid phase transition. This is observed in finite nuclei as the 
even-odd staggering of masses and is expected to occur in infinite nuclear matter inside neutron stars. Calculations 
t-H ■ in finite nuclei usually assume an effective pairing interaction fitted to the available data. Calculations in nuclear 
matter are using predominately the bare nucleon-nucleon interaction in the gap equation. 

Due to the short range repulsive core in the nuclear potentials Brueckner type ladder resummation of the interaction 
in medium is necessary On the other hand, the effective interaction in the gap equation should be two-particle 
irreducible, i.e. without resummation of the particle-particle (and hole-hole) ladder Q. Thus a good starting point 
could be the bare nucleon-nucleon interaction, the same as used in the ladder approximation for the two-particle 
f^) . correlation. 

The value of the pairing gap depends on nucleon single particle energies. Using single-particle energies obtained 
from Hartree-Fock or Brueckner-Hartree-Fock (BHF) calculation a different value of the paring gap is obtained. It 
■ can be understood as a modification of the effective mass, and consequently a change in the density of states at the 
Fermi energy. Obviously the value the effective mass is an inherent part of effective parameterizations of the pairing 
interactions used in calculations in finite systems. 

Only relatively few works discuss possible in medium modifications of the pairing interactions beyond a change in 
the effective mass. Polarization corrections to the bare nucleon-nucleon interactions were discussed in Refs. It was 
found that the screening reduces the pairing gap by a factor ~ 3. Another in medium modification was analyzed in the 
framework of self-consistent nuclear matter calculations [||. Self-consistent nuclear matter calculations use off-shell 
nucleon propagators in the resummation of ladder diagrams for the self-energy. In this way a self-consistent spectral 
function can be obtained. It was found that the use of full spectral functions in the gap equation leads to a strong 
reduction of the critical temperature and of the superfiuid gap in comparison to the quasi-particle approximation. 
It is the goal of the present paper to identify the main cause of this modification of superfiuid properties of nuclear 
matter and to propose a renormalization of the quasi-particle approximation for the gap equation. In Sect. II we 
discuss modifications of the two-particle correlations in the normal phase. This enables us to calculate the critical 
temperature, at which long range two-particle pairing correlations appear. The critical temperature is estimated using 



full spectral functions and in the renormalized quasi-particle approximation. In Sect. Ill we present a discussion of 



the gap equation with full spectral function and compare it to its quasi-particle limit and to the usual mean-field gap 



equation. In the concluding section [V we identify the most important renormalization of the pairing interaction due 



to off-shell propagation and indicate its consequences for realistic nuclear matter calculations. 
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II. QUASI-PARTICLE LIMIT OF SELF-CONSISTENT LADDER RESUMMATION 



Nuclear medium is a relatively dense system of particles strongly interacting on short distances. Brucckner resum- 
mation of particle-particle ladder diagrams defines the so called in medium G-matrix 

< P |G(p > n)|p'> = v(p,p')+ ( ^ V {j>,4 {l ~ f ^ m ~ fM) <q|g(P,n)|p > , (i) 

where pi,2 = P/2 ± q. G-matrix resummation allows to define single particle energies and gives relatively good 
results for the saturation properties of nuclear matter. In the above equation and in the following we skip the spin, 
isospin indices which are implicitly summed over. The above equation corresponds to a resummation of particle- 
particle ladders, with medium effects entering through the Pauli blocking factors 1 — f{u) p ) in the numerator and 
single-particle energies tu p in the denominator. Most advanced calculations in the Brueckner scheme use the so called 
continuous choice for the single particle energies u> p , self-consistently defined by the G-matrix m. 

Another approach starts from the T-matrix approximation for the two-particle correlations £fl . In this scheme the 
ladder diagrams include both particle particle and hole-hole propagation. The Pauli blocking factor (1 — f{u pi ))(l — 
f(tu P2 )) in the G-matrix equation is replaced by 1 — f(tu Pl ) — f{u) V2 ) in the equation for the retarded T-matrix 

<p|T(p,fi)| P '> = v(p,p')+ [ ^* V (p,k) {1 - fM < qir(p,n)| P ' > . (2) 

J (2ny i 1 l - u Pl — U! P2 + te 

The imaginary part of the retarded self-energy in the T-matrix approximation is 

f d 3 k f \ 

ImS(p,u) = J — ^ < (p-k)/2|/mr(|p + k|,w fc + w)|(p-k)/2>^ (f(w k ) + b(u + w k )) , (3) 

where b(u>) is the Bose distribution. < . . . >a denotes antisymmetrization of the T-matrix (also in the spin, isospin 
indices not explicitly show) . The real part of the self energy consists of the Hartree Fock self-energy and a dispersive 
contribution obtained from ImE 

ReH M =H HF{p)+ V j^JI^l. (4) 



The imaginary part of the self-energy is usually neglected leading to the quasi-particle approximation for the two- 
nucleon propagator in the T-matrix (Eq. |^). 

Allowing for off-shell propagation of nucleons and taking the self-energy self-consistently (also its imaginary part) 
requires the use of full spectral functions in the calculation resulting in a more complicated expressions for the T-matrix 
and the self-energy M 



<p|T(P,fi)|p > = V(p,p) 



1 ' ^(p,q)^— — -- L A(p 1 ,uj 1 )A(p2,uJ2) < q|T(P,fi)|p 5(5) 



and 



2ir J 2ir J (2ir) 3 ' fl — u>x — u>2 + ie 
du}\ f d 3 k 



ImE+(p,uj) = j ^ J ^A(k,Lu 1 )<(p-k)/2\ImT(p + k,Lu + u Jl )\(p-k)/2> A + g{u + .(6) 

Equations (ph, (pi) and (0) have to be solved iteratively with the constraint on the assumed density p p, M 



duj f d 3 p 
2^ J (2tt) 3 



p= I — I T^^Mm ■ (7) 



Results of the self-consistent calculation and of the quasi-particle approximation are very different. In medium cross 
sections for nucleon-nucleon scattering are smaller when using off-shell nucleons [0. Also the critical temperature is 
strongly reduced when using full spectral functions 0. 

The T-matrix approximation for the two-particle propagator is directly related to the superfluid gap properties by 
the Thouless criterion || and the condition for long range order jl(| ^l). On the other BHF calculation using the 
G-matrix are are much developed for realistic interactions. In this work we derive an improved gap equation which 
could use the results of advanced G-matrix calculations as an input. 
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A. Thouless criterion for superfluidity 

The critical temperature can be obtained from the Thouless criterion for superconductivity 01. This is the tem- 
perature where a singularity in the T-matrix appears at twice the Fermi energy (CI = 0) and zero total momentum 
of the pair (P = 0). It means that the real part of the inverse T-matrix develops a zero-eigenvalue at the critical 
temperature 

/ < Pl ReT_1 ( P = °> U = °)IP' > A (P) = ■ ( 8 ) 

It is equivalent to the existence of a nontrivial solution of the gap equation at T c 
A( P ) = 

d 3 k fduO fdu\ v , A(k,L,-Lj')A(k,L,)(l-f(LJ-Lj')-f(Lj')) 



(2tt) 3 J 2tt J 2tt v ' u 

The two propagators in the gap equation above enter with the full spectral function. On the other hand the BCS 
quasi-particle gap equation is 



(2nf K(P ' k) 2( k 



A (P) + / V(P, V- K? A ( k ) = (10) 



with £ p = p 2 /2m + Yihf{p) — P, and it corresponds to the Thouless criterion for the quasi-particle T-matrix, i.e. the 
appearance of a singularity in the T-matrix given by Eq. (^) . Because of additional averaging over spectral functions 
in the gap equation with off-shell propagators, a different value of critical temperature comes out. 

The self-consistent T-matrix calculation was done for a simple, S-wave, interaction and compared to the quasi- 
particle approximation for the gap equation ||. At a density of 0.45 of normal nuclear density po it was found that 
the critical temperature was reduced from T c = 5MeV in the mean-held gap equation with Hartree-Fock single particle 
energies to T c = 1.6MeV in the T-matrix approximation with off-shell propagation ||. Below T c a modified T-matrix 
resummation was used in Ref. ||. The resulting superfLuid gap is significantly smaller than the one obtained from 
the usual gap equation with quasiparticles. Modifications to the gap equation below T c coming from the use of full 



spectral functions will be discussed in Sect. III. Here we concentrate on the T-matrix equations in the normal phase, 
which is sufficient for the calculation of the critical temperature. The simple Yamaguchi S wave interaction Q does 
not permit calculations at low temperature for normal density nuclear matter. 

B. Renormalized quasi-particle interactions 

Excitations in the Fermi liquid close to the Fermi energy can be described by quasi-particles. Quasi-particles are 
propagating on shell with dispersion relation modified by the presence of the medium. Also the scattering amplitudes 
between quasi-particles are modified by the medium. One of these modification comes from the quasi-particle limit 
in the propagator of two nucleons jl], Q . Two particle Green's function can be formally written as a resummation 
of particle-particle (and hole- hole) ladder diagrams starting from in medium two-particle irreducible vertex jl], [l3) . 
For short range interactions the two-particle irreducible vertex can be approximated in the lowest order by the 
bare nucleon-nucleon interaction, leading to the T-matrix equation (|^). It should be pointed out that due to the 
instantaneous form of the interaction in the T-matrix equation, the full T-matrix depends only on the total energy, 
and its equation takes a simple form for the retarded T-matrix (Eq. [5]), also at finite temperature. 

Close to the Fermi energy the spectral function becomes peaked around the quasi-particle pole 

A(p,w) = Z p 2t:S(lo - co p ) + R(p : uj) , (II) 

where 

Ml (12) 

and R(p,uj) is the regular part, smooth in the vicinity of the quasi-particle pole. The retarded propagator of two 
nucleons appearing in the T-matrix ladder can be written as 

(l - /(wi) - /(w 2 )) 

±— — J -A(p 1 ,co 1 )A(p 2 ,co2) = B re9 {p ll uj ll p 2 ,Lo 2 ,d) 

\l — LO l — U>2 + 
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+ (2tt) 5(0;! -o; pi )d(a; 2 -u; p J ^- - . (13) 

_B re9 describes the part of the propagator of two nucleons which cannot be written as propagation of two quasi-particles, 
it originates from the background part R(p, ui) of the spectral function. Besides this contribution, the propagator of 
two-nucleons differs from the one used in the quasi-particle T-matrix (^) by the presence of renormalization factors 
Z p (Eq. |l2|) . The T-matrix with quasi-particle propagators takes the following renormalized form 

<p|T(p,fi)| P ' > = < P \v ren (p,rt)\p > 

+ fwi> < pr en (*^)lq > ^^ (1 Q ' /Kl) ' /( ? 2)) < > , (w) 

J (^r u-u Pl - u P2 + it 

with the renormalized interaction y ren given by 

<p|V r " l (P,n)|p >=V(p,p') + J^-J^J ||%q)B(pi, W i,P2^2^) < q\V ren (P,n)\p > (15) 

The bare interaction is renormalized by contributions from background parts of the spectral functions, this involves 
integration over energies far from the quasi-particle pole. The background part of the spectral function is of course 
necessary to recover sum rules for the particle strength. However, the full treatment of the renormalized interaction 
is difficult and in the following we replace it by the bare interaction y ren ~ V in the gap equation. The only remnant 
of the dressing of nucleons in medium are the single particle energies 

Up = £p + ReT,(p, lu p ) (16) 

and the Z p factors in the homogeneous term in the T-matrix equation. It will turn out that these are the dominant 
modifications responsible for shifting the critical temperature. In the vicinity of the pole the T-matrix equation is 
dominated by the homogeneous term. In this region the interaction for quasi-particles is effectively renormalized by 
a factor Zp^ Zp 2 . 



C. Two-particle pole with renormalized interactions 

According to the Thoul ess c riterion, a pole in the T-matrix at the Fermi energy (CI — 0) means that superfluid long 



range order sets in (Sect. [I A). As mentioned earlier, standard mean-field gap equation uses the following kernel 



(l-2/(Cfe)) 

V(p,k)± >- (17) 



2Cfe 

with mean-field single particle energies Ck = k 2 /2m + Y^HF(k) — H- It leads to a critical temperature T c — 5MeV. 
One could take single particle energies ui p beyond Hartree-Fock, e.g. from BHF or T-matrix calculations. However, 
the difference between the resulting effective masses is of the order Z p . Hence we stay at the order of Hartree-Fock 
effective mass in the standard gap equation, including the modified single particle energies ( |l6| ) only together with 
Z p factor renormalization. 

The second estimate for the critical temperature is obtained from the condition of appearance of the pole in the 
T-matrix with full self-consistent spectral functions. It is equivalent to the following kernel in the gap equation 

W^ y(p ' k) « ' (18) 

where A(p, u>) is obtained from full self-consistent calculation of normal nuclear matter at finite temperature with 
off-shell propagators in the T-matrix ladder ||. It gives a very different value for T c = 1.6MeV at p — A5po- 

Finally we can use the renormalized interaction strength in the homogeneous term of the T-matrix equation (Eq. 
[l4]) . It is equivalent to using a renormalized interaction in the kernel of the quasi-particle gap equation 



1 - 2/M) 



y fr'*>^ V 2u k ^ (19) 

with the single particle energies lu p and Z p factors obtained from the full self-consistent T-matrix calculation (Eq. 
|l6|) . It gives a value of T c = 2.2MeV, much closer to the result of the calculation with full spectral functions. The 
Z factor obtained from the full self-energy is not close to 1, at the Fermi energy we find in our model calculation 
Z p , ~ 0.7. Such a small value of the renormalization factor can explain the large difference in critical temperatures 
found in the standard quasi-particle gap equation and in the one with full spectral function. 
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III. GAP EQUATION WITH DRESSED PROPAGATORS 
A. Gap equation with full spectral function 

Below the critical temperature a nonzero solution of the gap equation is possible. The kernel of the gap equation 
is very similar to the two-particle propagator in the T-matrix resummation (Eq. ^|) . However one of the propagators 
includes the normal self-energy as well as the off-diagonal one A(p) ||. The full retarded propagator can be expressed 
using the normal propagator G(p, lo) (which includes only the normal self-energy) 

G ^ ) = GM-^ + \m P) G^,-.)- (20) 

As a result in the kernel of the gap equation 

A( P ) = - [** flOL [ ^y ip , k) A ^~ ^ A ^)(i - f± - J) - /(„')) 

J [2tiY J 2tt J 2tt lo 

two spectral functions A(p, lo) = —2ImG(p,Lo) and A s (p, lo) = —2ImG s (p,io) appear. Both spectral functions can be 
calculated knowing the two self-energies £(p, lo) and A(p). The off-diagonal self-energy, ie. the superfluid gap can 
be obtained from the gap equation ( pl| ) . On the other hand the normal self-energy cannot be calculated in the usual 
T-matrix approximation as above T c . This is related to the appearance of the Cooper instability in the two-particle 
propagator in the ladder approximation. 

The T-matrix equation can be modified by introducing also anomalous propagators in the ladder pi] . In Ref. || 
we used a simpler modification of the ladder diagrams, using a mixed ladder with one full propagator and one with 
only normal self-energy included. This means that the singularity of the T-matrix which appears at T c at the Fermi 
energy O = and zero total momentum of the pair, stays there also below T c . This reflects the presence of long 
range order |ll). In such a way we were able to obtain self-consistent solutions for the normal self-energy in the 
ladder approximation and for the superfluid gap in the mean-field approximation, using at all stages full spectral 
functions A(p, lo) and A s (p, lo), without quasi-particle approximation. The results for the superfluid gap at several 
temperatures below T c are indicated by triangles in Fig. [l. This calculation cannot be extended straightforwardly to 
zero temperature because at the Fermi energy the imaginary part of the self-energy vanishes and the discretization of 
the spectral functions is not possible close to the Fermi energy. At the Fermi momentum quasi-particles appear with 
small width. 

By comparing the imaginary part of the self-energy at two temperatures T = 1.7MeV and 1.2MeV one notices 
that the only modification is the reduction of single-particle width close to the Fermi energy with temperature, in 
agreement with general properties of Fermi liquids. Thus we can take as an approximation a temperature independent 
imaginary part of the self-energy. We proceed by taking the imaginary part of the self-energy ImT,(p, lo) as calculated 
at T = 1.63MeV, slightly above T c . Obviously the dispersive contribution to the real part of the self-energy is also 
fixed (0). The Hartree-Fock energy is obtained from 



/d 3 k 
_V(|p_k|/2,|p-k|/2)n(p) (22) 



where the momentum distribution is given by 



n(p) = J ^-A s (p,Lo)f(Lo) . (23) 
The shift of the Fermi energy from the value lo = at T — 1.63MeV to keep the density 

p = - 45po = / Wf n{p) (24) 

constant when decreasing the temperature to zero is only 0.2MeV and has negligible effect on the Hartree-Fock energy. 
Thus the main change in the self-energy when decreasing the temperature occurs in the off-diagonal part A(p). At 
any given temperature the superfluid gap is obtained from the gap equation with full spectral function (|2l|). Where 
the dependence of the kernel of the gap equation on A(p) enters through the spectral function 

A s (p,lo) = -2 ( (lo + C P + ReS+(p, -w)) 2 Im£+(p, lo) + ImE+(p, -lo)A 2 { P ) + (lm£+(p, -Lo)flmT,(p, lo) ) / 
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FIG. 1: Superfluid gap at the Fermi momentum as function of the temperature at p = A5po- Points denote results obtained 
from self-consistent solution of the ladder self-energy and gap equations. The solid line represents results obtained using the 
imaginary part of the self-energy fixed at T = 1.63MeV and solving self-consistently the Hartree-Fock and the full gap equation 
(21). The long dashed line is the mean- field quasi-particle solution of the gap equation with Hartree-Fock single-particle energies 
(26). The dashed line is the result of the quasi-particle approximation (|6|) without renormalizing the superfluid gap (jjif). The 
dashed-dotted and dotted lines are the results of the quasi-particle approximation (^6|) with renormalization of the superfluid 
gap by Z p (^) and by 1/(1 - 0(p,u p )) (^) respectively. 



- £ P - Re£+(p, w))(w + t P + Re£+(p, -u)) - ImS+(p, w)Im£+(p, -w) - A 2 (p)j 2 
Im£ + (p,w)(w + £p + ReE + (p,-a;)) + Im£+(p, - £ p - Re£+(p, a;))) ] . (25) 



In the above equation TtoE is fixed and temperature independent, according to our approximation, and i?eS depends 
only very weakly on the temperature and the superfluid gap through the Hartree-Fock energy. The solution of the 
gap equation using a fixed single particle width down to zero temperature is represented by the solid line in Fig. [I]. 
As expected close to the temperature where the imaginary part of the self-energy was fixed (T = 1.63MeV) the result 
is close to the fully self-consistent solution denoted by the triangles on the figure. As the temperature is lowered 
some deviations appear. It is due to the decrease of the single particle width close to the Fermi energy in the fully 
self-consistent solution (Fig. ||). In the following we will compare the solution of the gap equation with full spectral 
function to the corresponding quasi-particle approximation. Since the imaginary part of the self-energy is fixed at 
T = 1.63MeV, the properties of the quasi-particle pole (oj p and Z p ) are taken from the self-consistent solution at the 
same temperature. Our aim is to explain the big difference between the mean-field solution with on-shell propagators 
(long-dashed line in Fig. [j]) and the gap equation with full spectral function below T c (solid line). For this purpose it 
is sufficient to compare the approximation with nontrivial but fixed spectral properties with its quasi-particle limit. 
The mean-field BCS gap equation is 



d 3 k n-2f(E k 
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FIG. 2: Imaginary part of the self-energy at zero momentum for two temperatures T = 1.7 (solid line)and 1.2 MeV (dashed 
line) . 

where Ef. = -^/Cf + A(fc) 2 . It gives significantly larger values for the superfluid gap (long-dashed line in Fig. [j]) than 
the full solution (solid line) in all ranges of temperatures below T c . 



B. Quasi-particle approximation for the superfluid spectral function 

Above T c we have used the standard form of the quasi-particle approximation for the spectral function A{p 1 uj) — 
2u Z p 5{uj — uj p ). Using it we have found a significant reduction of the critical temperature with respect to the mean-field 
approximation. To obtain a quasi-particle approximation for the kernel of the gap equation ( |2l| ) with full spectral 
function we have to construct an approximation for the superfluid spectral function A s . We can write the spectral 
function ( p5|) putting an infinitesimally small imaginary part of the self-energy (ImT,(p, uj) = —irf) 

A s (p, uj) = -11m ^ (uj + £ P + ReE(p, -uj) + ir^ / ({uj - £ v - ReE{p, uj)) 

(w + ^ + i?eS(p,-w)) + A 2 (p)+i77(2w-i?eE(p,w) + i?eI](p,-w)))^ , (27) 

with £ p = p 2 /2m — /i. It is useful to define the even S(p,uj) = i?eS(p, uj) + i?eS(p, — uj) and odd luO(j>, uj) = 
Re£{p,uj) — i?eS(p, —a;) parts of the real part of the self-energy with respect to the Fermi energy |T^]. Notice that 
we can take S(p, 0) = redefining fi. We have 



A s (p,uj) = -2Im\ 



uj{l-0{p,uj))+i P + S{p,uj) 



((1 - 0(p, uj)) - E{p, uj)) {uj{1 - 0( P , uj)) + E(p, uj)) + i V w{l - 0(p, uj)) 



(28) 



where 



E( P ,uj) = tJ(Z p + S(p,uj)) 2 + A 2 (p) . 
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The above expression has poles at 

E(p, e p ) E p 



p l-0(p,e p ) l-O p ' 
which gives two quasi-particle contributions to the spectral function on both sides of the Fermi energy 



A s (p^) = 2^z'l Ep+ 2 ^ +Sp S(u;-e p ) + Ep ^ Sp ' 5(u> - e p )\ , (29) 



with S p = S(p, e p ) and 

/ d(E(p,u)/(l-0(p,u)) 
K 1 = f 1 - ^ -|«=%, I (1 - O p ) . (30) 

It will be useful to relate this new pole renormalization strength to the usual renormalization factor Z p (Eq. |l^) . We 
can write 

/ _ x _ , _ \ j E-/ p sign(p - Pf )(£ p + S p ) dSjp, w) 

with sign(a;) = Q(x) — ©(— x) . For small values of the superfluid gap the superfluid quasiparticle position sign(p— Pf)tp 
is very close to the position of the pole in the normal propagator lo v = £ p + i?e£(p, u p ). Substituting lo p for the energy 
argument in the first term in Eq. (|3l]) we get 

7 ' -i 7-1 _ E p - signjp - p f )(Zp + S p ) dS(p,uj) f 

z p ~ z v ~ Ep d^ lu=e * ■ {32) 

The first factor on the right hand side of the above equation is small except close to the Fermi energy where the 
second factor precisely vanishes. Thus we expect that the new quasi-particle pole renormalization factor Z p can be 
approximated by the renormalization factor of the normal spectral function Z p corresponding to the same momentum. 
A numerical calculation confirms this to a very good accuracy (Fig. ^|). 
The spectral function takes the form 

A s (p, (J) — 2irZ p ^uIS(lo - e p ) + Vp5(uj - e p )^ , (33) 

with coherence factors 

U>1) = *,+(-) fe+MMp) . (34) 



v 



These coherence factors can be very well approximated by an expression similar to the one used in the usual mean-field 
gap equation 

ulivl) = £p + 2£ ( ; H ■ (35) 

The approximation works to within 1% for A(p/) < lOMeV. 

Substituting the quasi-particle expressions for the spectral functions A and A s into the kernel of the gap equation 
one obtains 

(36) 

We have obtained an expression for the gap equation with a kernel very similar to the mean-field BCS one (Eq. fl7| ), 
but with the interaction renormalized by Zp and a different quasi-particle energy e p . 

In the limit of vanishing gap A(p) — ► all quasi-particle approximations for the gap equation reduce to the 



Thouless condition for the renormalized T-matrix equation (Eq. 19). Accordingly we have the same condition for T ( 



as discussed in Sect. [I C 
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FIG. 3: Quasi-particle pole renormalization factor for the full spectral function Z v (Eq. p0[) as function of momentum (solid line) 
compared to the quasi-particle renormalization in the normal spectral function Z v (EqTTld) (dashed line) for A(pf) = 2MeV. 



C. Quasi-particle energies 



In order to relate our quasi-particle gap equation (pq ) to usual approaches we have to obtain an approximation for 
the quasi-particle energies 



Vfe, + S P ) 2 + A 2 (p) 



(37) 



To extract the positions ±e p of the poles of the spectral function one has to know the real part of the self energy for 
energies on both sides of the Fermi energy. We would like to obtain an expression which could be used as a correction 
to calculations using quasi-particle approximation like BHF. 

In BHF approaches one calculates a single particle energy which includes dispersive corrections to the position of 
the quasi-particle pole (Eq. |l6| ) in the usual spectral function A(p,oj). A simple approximation would be to use u> p 
instead of ( p in the expression for the quasi-particle energies in the super-fluid 



>2 + A*(p) 



(38) 



Clearly the above expression does not give the right value for the energy gap at the Fermi momentum (Fig. ^). The 
energy gap at pf should be 



A(pf) 7 w v 
e P p = 1 ~ n = z pf A \P f i ■ 



(39) 



Let us define a renormalized energy gap 



(40) 
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FIG. 4: Position of the superfiuid quasi-particle pole e p (solid line) as function of quasi-particle energy compared to the naive 
expression + A 2 {p) (dashed line) (Eq. ||) for A(p/) = 3MeV. The quasi-particle energies obtained using renormalized 

superfiuid gap y/ ujp + A(p) 2 /(1 — O p ) 2 (Eq. ptT) cannot be distinguished from the solid line on the scale of the figure. 

Results for the quasi-particle poles using this renormalization of the superfiuid order parameter 

e p = Juj + frip) (41) 

are indistinguishable from the solution ( |37| ) in Fig. [|. It should be noted that for the calculation of the factor (1 — O p ) 
we can take O p — 0(p,e p ) ~ 0(p,u p ). So that the approximation for e p can be expressed using the quasi-particle 
pole energy uj p and the values of the real part of the self energy at ±w p . However, usually in BHF calculations one 
does not calculate the factor 1 — O p , since the value and the derivative of the real part of the self-energy is known 
only at lo p and not at — uj p . In this case the position of the poles in the superfiuid can be approximated by 

e p = y/uj + A*(p) (42) 

with the renormalization of the superfiuid gap taken as 

A(p) = A(p)Z p . (43) 

A(p) is very close to A(p) around the Fermi energy. Some differences appear only around p — (for the assumed 
interaction) . but this region is less important for the solution of the gap equation. Again the energies obtained using 
expression ( [42| ) are indistinguishable from the solid line in Fig. ^. 

D. Results for superfiuid gap in quasi-particle approximation 



In previous sections we derived a quasi-particle approximation for the kernel of the superfiuid gap equation (Eq. 
11) . The results of the full gap equation with integration of over the energies in the fermion propagators (solid line 
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TABLE I: Values for the critical temperature and superfluid gap at zero temperature for different approximations. 



appi UA.imaiioii 


± c \ i vie v j 


^\Pf)t— l^ivie v ^ 


aCtii-.i it 


off-shell gap 


1.6 


4.3 


2.7 


equation (|2l|) 








quasi-particle app. with 


2.2 


6.0 


2.7 


A(p) (Eqs. pq, (41|) 








quasi-particle app. with 


2.2 


5.7 


2.6 


A(p) (Eqs. p6|, |42|) 








quasi-particle app. with 


2.2 


4.0 


1.8 


A(p) (Eqs. pq, p8p 








mean-field BCS 


5.0 


10.8 


2.2 


(Eq. pel) 








weak coupling BCS 


1.76 


weak coupling BCS with 






2.5 


renormalized energy gap 









in Fig. can be compared to quasi-particle results. The dotted line in Fig. |l| represents the solution of the gap 
equation ( |36| ) with renormalization of the energy gap A(p) for the calculation of the superfluid quasi-particle energies 
(Eq. [II]). The dashed-dotted line represents an analogous calculation but with the energy gap A(p) (Eq. ^). Both 
quasi-particle approximation are very close to each other with critical temperature T c — 2.2MeV. The value of T c is 



the same as obtained in Sect. [I C using renormalized interactions in the ladder approximation. This value of critical 
temperature is much closer to the one obtained from the self-consistent solution T c — 1.6MeV than was the mean- field 
value T c = 5MeV. Also the value of the superfluid gap is much closer to the one calculated using full of-shell spectral 
functions (Table |jj). The superfluid gap at zero temperature in the quasi-particle approximation A(pp) — 6MeV is 
much closer to the solution of the off-shell gap equation A(pp) — 4.3MeV than the one obtained from the mean-field 
BCS solution A(pp) = 10.8MeV. From the results in Table | we can notice that the ratio between the superfluid 
gap at zero temperature and the critical temperature is larger in the solution of the full gap equation ( ^l|) and in 
quasi-particle approximations ([56]) with renormalization of the energy gap ( f4l| or |4^) , than in the mean-field solution 
and in quasi-particle approximations without renormalization of the superfluid energy gap ( |36| , |38| ) . In the weak 
upling limit we have |l4|, [l5j 

A(P ; )|T - = ne^ ~ 1.76 . (44) 



If the energy gap A(jpp) is renormalized with respect to the superfluid gap A(pp) (the off diagonal self-energy) we 
have 



A(p F )|r=o = A(pf)| r = (l-Op F ) = A( PF )\ T=0 ^ L76 

-M3 T-C 2c F PF 



2.5 . (45) 



Ratios of the superfluid gap and of the critical temperature for the solution of the full gap equation and quasi-particle 
approximations with renormalized superfluid gap is close to 1.76 /Z PF (Table |). On the other hand mean- field gap 
equation and quasi-particle approximation without renormalization of A(p) give A{pf)\t=o/T c closer to 1.76. These 
relation can be fulfilled only approximately since we are far from the region of applicability of the weak coupling BCS 
solution. 



IV. CONCLUSIONS 



Numerical solution of the gap equation with full spectral functions ( pT| ) shows a strong reduction of the superfluid 
gap and of the critical temperature with respect to the mean- field BCS solution (^6|). To understand this effect we 
constructed a quasi-particle approximation for the full gap equation. The effects of nontrivial spectral functions can 
be approximated using a renormalized strength of the interaction V(p, k)Z^. Also the energy gap in the calculation of 
the quasi-particle poles in the superfluid is renormalized. For the renormalization of the gap we used two expression 
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A(p)/(1 — Op) and A{p)Z p . In both cases we obtained results for A(p) and T c much closer to the full solution than the 
mean-field approximation. Having at one's disposal only the single particle energies uj p and the Z p factors obtained 
from a realistic BHF type calculation, the gap equation can be corrected. A reduced interaction strength V(p, k)Z% 
must be taken and a factor Z p appears between the energy gap A(p) and the off-diagonal self-energy A(p). This 
is equivalent to solving a gap equation for A(p) with reduced interaction V(p, k)Z p Zk- In the weak coupling limit 
A(p)| T=0 ^ 1.76T C but A(p)| T=0 ~ 1.76T C /Z PF . 

In the illustrative model here presented the scattering corrections are very strong (Z PF ~ .7) and there are still 
some differences between the improved quasi-particle approximation and the full solution of the gap equation. We 
expect that at normal nuclear density where the Z factor is closer to 1, quasi-particle approximation with renormalized 
interaction strength would be much closer to the full solution. This is the case for neutron matter where effects of 
renormalization of quasi-particle poles and shifts in single-particle energies are smaller Jl6| . In future work we plan 
to investigate the effects of the renormalization of the interaction by background corrections (|l5|) . 

In this investigation we used the bare potential for the two-particle irreducible vertex. However, it is known that 
polarization effects reduce the superfluid gap by a factor ~ 3 (2j . We must conclude that although there are no ladder 
corrections to the interaction in the gap equation, other many-body effects modify the effective interaction. Due to 
the exponential dependence of the gap solution on the strength of the interaction, these usually neglected corrections 
modify strongly superfluid parameters in nuclear matter. 
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